Linearized Tensor Renormalization Group Algorithm for Thermodynamics of Quantum Lattice 

Models 
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A linearized tensor renormalization group (LTRG) algorithm is developed to calculate the thermodynamic 
properties of low-dimensional quantum lattice models. This new approach employs the infinite time-evolving 
block decimation technique, and allows for treating directly the transfer-matrix tensor network that makes it 
more scalable. To illustrate the performance, the thermodynamic quantities of the quantum XY spin chain as 
well as the Heisenberg antiferromagnet on a honeycomb lattice are calculated by the LTRG method, showing 
the pronounced precision and high efficiency. 
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Since the appearance of White's density-matrix renormal- 
ization group (DMRG) theory [1], the numerical renormaliza- 
tion group (RG) approaches have achieved great success in 
studying low-dimensional strongly correlated lattice models 
In the past few years, a number of RG-based methods, 
ej*., the coarse-graining tensor renormalization group (TRG) 
Of0|, projected entangled pair states jfj, entanglement renor- 
malization |7], the infinite time-evolving block decimation 
(iTEBD) dsnT finite-temperature DMRG |S O, etc., have 
been proposed inspired by the quantum information theory. 
In spite of the great success in one- and two-dimensional (ID 
and 2D) lattice models, it is still quite necessary to develop 
new algorithms to improve the accuracy and efficiency of nu- 
merical calculations for strongly correlated systems. 



In this Letter, we propose a new algorithm to simulate the 
thermodynamics of low-dimensional quantum lattice models. 
Our strategy is first to transform the D dimensional quantum 
lattice model to a D + 1 dimensional classical tensor network 
by means of the Trotter-Suzuki decomposition Hill , and then 
to decimate linearly the tensors following the lines developed 
in the iTEBD scheme to obtain the thermodynamics of the 
original quantum many-body system. This algorithm is so 
dubbed as the linearized TRG (LTRG). As is known, the pre- 
vious real space TRG approach deals with the 2D tensor net- 
work with exponential decimation in the coarse-graining pro- 
cedure, which was shown effective for both 2D classical and 
quantum lattice models S [l2T - [l5ll . For the best illustra- 
tion of the algorithm and performance of the LTRG approach, 
we take the exactly solvable ID quantum XY spin chain as a 
prototype. The results show that the precision of the LTRG 
method is comparable with that of the transfer-matrix renor- 
malization group (TMRG) [ 16], the method that is quite pow- 
erful for simulating the ID quantum lattice models at finite 
temperatures (e.g. Refs. Jl3,H)- To demonstrate its scalabil- 
ity, a LTRG result with remarkable precision for a 2D Heisen- 
berg antiferromagnet on a honeycomb lattice is also included. 



Let us start with the Hamiltonian of a ID quantum many- 




FIG. 1: (Color online) (a) A transfer-matrix tensor network, where 
each bond denotes the a index in Eqs. (O and l[3}. (b) A local 
transformation of a fourth-order tensor into two third-order tensors 
through a singular value decomposition (SVD). (c) Transform the 
transfer-matrix tensor network to a hexagonal one. (d) By contract- 
ing the intermediate bonds marked by dashed ovals in (c), one gets a 
brick wall structure with the 4th-order tensors in the bottom line. 



body model given by 

A' 

H = J]h u+l =H l+ H 2 , 
;=i 

N/2 Nil 

H\ = ^jhzi-\,2i, H 2 = ^%2!+i, (1) 



where N (even) is the number of sites. By inserting 2K (large 
K) complete sets of states {\o J i )}(o J j = 1, ■ • • ,D) with i the 
site index and j the Trotter index, the partition function of this 
model can be represented as 
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FIG. 2: (Color online) A local evolution of the tensors by contraction 
and SVD. (a) Contract the intermediate bonds; (b) obtain a 6th-order 
tensor 0; (c) calculate the singular value decomposition (SVD) of 
O, and update the tensors M a h and A. The above manipulation has a 
computational cost that scales as 0(D 6 Dl). 



Z N - Tr [e- m,K e-^ /K ] K 

x (o»...o*j\eW'\a>»\.Sf 1 ), (2) 

where the periodic boundary conditions along both spatial 
and temporal directions are assumed, i.e., cr\ = cr 2K+l and 
cr' x - o~ J N+l . Since the terms within H\ (and Hj) mutually com- 
mute, Eq. (|2]i can be further decomposed as 

N/2 K 

Zn * Z 11 II v <^vrvM v 4°-iv°-T i °%:i • (3) 

(o-f) «=1 J= l 

where the transfer matrix, v cr ^ iyCr ^ cri = 
{a\(Tt\evp(-lihij+\[K)\(r20-}), is a 4th-order tensor. Ob- 
viously, the partition function, Eq.©, can be viewed as a 
classical transfer-matrix tensor network, as illustrated in Fig. 
1(a). 

The partition function can be obtained by summing over 
all the intermediate states \cr J ), namely, contracting all the 
bonds <x in the tensor network. This procedure is accom- 
plished by first making a singular value decomposition (SVD) 
of v-tensors in the following way 

D 2 

x=l 
D 1 

= ^^J^ a ) xo - x0 - 2 {Tjy)x,cT^,cr4t (4) 

.v=l 

where the diagonal matrix A collects D 2 singular values, 
and two auxiliary tensors (T a ) Xt0 - u(T2 = U a - xa - 1 ^yfA~ x and 
(Tb)x,o-,,o; = Va- 30 - 4tX are introduced for convenience. Af- 
ter this transformation, the square tensor network becomes a 
hexagonal one with two 3rd-order tensors T a and Tb, as de- 
picted in Fig. [TJb). Then, one contracts the cr-bonds encircled 
by the dashed oval lines between the last two rows in Fig. [TJ C X 
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FIG. 3: (Color online) An successive projection of each row of ten- 
sors onto the MPO in the bottom line [(a)-(c)]. After the projection 
along the Trotter direction, by tracing out the physical indices t and 
b of the MPO, one may get a ID matrix product, of which the trace 
can be obtained by a matrix RG procedure [(d)-(g)]. 

which leads to the two 4th-order tensors 

D 

v=l 
D 

(Mb)/3,t 2 ,y,b2 - '*^(T a )y^ tl {T b )f3^ bl , (5) 

z=i 

which form a matrix product operator (MPO) lying in the bot- 
tom line of the whole tensor network, that can also be viewed 
as a "superket" in the operator Hilbert space [19]. Each hor- 
izontal bond between M a and Mb is assigned with a diagonal 
matrix A\£. Finally, we obtain a tensor network with brick 
wall structure as shown in Fig. [T](d). 

Next, one can project the tensors T a ,b onto M a ,b succes- 
sively. At each time, we project one row of tensors T a and 
Tb followed by updating M a ,b and A\ t 2- After two projections, 
the system evolves one Trotter step forward. This procedure is 
illustrated in Fig. [2] One first contracts the cr-bonds between 
M-tensors and T-tensors in Fig. [2] (a) to obtain a 6th-order 
tensor in Fig. |2](b) 

Oy, a ,bi,z,y,b 2 = / \ (^l)g {M-a)aMfiJbi (Affc)jS,t 2 ,y,& 2 (^1)7 

(T a ) XJl ,y(T b ) x , z , l2 , (6) 

and then, takes a SVD of the (9-tensors (after matriciza- 
tion). O ya b u zyb 2 - UyahfiitypVj^, while keeps only 
the largest D c singular values of A' T One can define new 
M-tensors (M' a ) ajJ3M = Uy ablJ3 /(Ai) a and (M' b )p A y M = 
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V zy b 2 fi/(Ai) y , and update the horizontal bonds with A' T Af- 
ter these operations, the last row of the tensor network is half 
updated as shown in Fig. |2](c). To project the next row of ten- 
sors, one can simply exchange M a and Mj as well as A\ and 
/I2 in Eq. ©. These two successive projections make up of a 
full Trotter step r, as illustrated from Fig. [3ja) to[3fc). In each 
Trotter step, the transfer-matrix tensor network is decimated 
linearly with only 0(D C ) singular values discarded, which im- 
proves greatly the efficiency compared with the original TRG 
approach where 0(D") (n — 2 for honeycomb network) ones 
are discarded in the coarse-graining procedure M20I1 . 

In order to avoid the divergence in the imaginary time evo- 
lution, one has to normalize all the singular values in A with 
its largest one n, in z'-th step. After projecting all the 7-tensors 
at inverse temperature /3, one is left with the matrix product 
density operator of the present system. It consists of 4th-order 
M-tensors [see Fig. [3(c)], each of which has two legs with 
physical indices t and b in the Trotter direction, that can be fur- 
ther traced out due to the periodic boundary condition. Thus, 
we obtain a ID matrix product (MP) extended in the spatial 
direction, where the matrices are labeled as cM a ^ as shown 
in Fig. [3](d). It is convenient to assume the number of ma- 
trices is 2 P . To get the trace of the product of these 2 P ma- 
trices, one can contract the neighboring matrices pairwise to 
obtain a new product of 2 p ~ l matrices, each of which should 
be normalized by the absolute value of its largest elements to 
avoid divergence. This contraction procedure is represented 
in Figs. [3jd)-[3jg). After p steps, the 2 P matrices shrink to a 
single one, of which the trace can be easily calculated. In each 
coarse graining step, all the normalization factors denoted by 
171 j with j = 1 , • • • , p need to be collected for the calculation 
of physical quantities, e.g., the free energy per site / at inverse 
temperature f3 = Kt can then be determined by the normaliza- 
tion factors n/s and m/s 

. 2K-2 p 

P 4=1 j=i 

. 2K-2 . p , 

!=1 j=\ 

In the above descriptions, we illustrate the LTRG algorithm 
by first decimating the tensors along the Trotter direction, and 
then contracting the matrices in the spatial direction. Alterna- 
tively, one can also perform the decimation first in the spatial 
direction, and then do the matrix contraction in the Trotter di- 
rection. 

As an example, we are going to demonstrate the efficiency 
of the LTRG algorithm by computing the free energy and 
other thermodynamic quantities of the quantum XY spin- 1/2 
chain with a local Hamiltonian ha+\ = -J(S x S j: l , + S y S y , ,) 
in Eq. (Q} with 7 = 1. We take the chain length to be 2 100 , 
which definitely reaches the thermodynamic limit. 

In Fig. HI we show the relative error of the free energy / 
with respect to the exact solution, i.e., Sf = \{f — f exact) I f exact], 
for different Trotter steps r = 0.1,0.05,0.02,0.01. We observe 




FIG. 4: (Color online) The relative error of the free energy per site, 
Sf, of the quantum XY spin chain at high temperatures. Sf converges 
rapidly with D c , and the lines with D c = 100 and 150 coincide with 
each other (t = 0.05, 0.02). In addition, the TMRG results (t = 
0.1, 0.05) are also presented for a comparison. 



that the accuracy is enhanced with decreasing t, as well as in- 
creasing D c . Owing to the close relation between iTEBD and 
DMRG, the truncation parameter D c plays a role similar to 
the number of states kept M in the TMRG method. As shown 
in Fig. H] we compare the LTRG results to those of TMRG, 
both of which show the same accuracy for r = 0.1 and 0.05. 
It is also noticed that the relative errors saturate rapidly with 
increasing D c , implying that the errors at high temperatures 
{e.g. T > 0.27) mainly originate from the Trotter-Suzuki de- 
composition. In order to check the truncation error, the LTRG 
algorithm is also tested at very low temperatures. In Fig. [5] 
the temperature is down to I = 7/120 with a Trotter step 
t = 0.05. As shown in Fig. 0(a), the accuracy of low T 
results is remarkably improved by increasing D c , and the rel- 
ative error 5f 7 x 10~ 6 at f3 = 120 for D c = 150. 

Besides the free energy, other thermodynamic quantities, 
such as the internal energy, can also be obtained. There are 
at least two ways to get them, one can either introduce some 
impurity tensors in the tensor network (see, for instance, Ref. 
8), or do a numerical differentiation of free energy with re- 
spect to temperature. Both ways are found to have a similar 
accuracy. In Fig. |5](b), the energy per site, e, is presented. 
We apply the LTRG algorithm to approach the ground state 
energy eo, and find the difference (e - eo)/<?o is about 10~ 4 at 
/? = 120 for D c =150, suggesting that the LTRG result is very 
close to the exact solution. The TMRG results with various 
M (up to M — 200) are also included in Fig. |5]for a compari- 
son. The relative errors for the free energy and internal energy 
are found to be of the same order down to ft - 120 for both 
approaches. 

The specific heat of the quantum XY spin chain is also 
calculated, as shown in Fig. |6ja). The LTRG results agree 
quite well with the exact solution both at high and low tem- 
peratures. As indicated in the inset, the accuracy will be en- 
hanced by increasing D c . For D c - 150, the LTRG results 
coincide with the exact solution down to very low tempera- 
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and boson models, and also provides a quite promising way to 
simulate the 2D quantum lattice models without involving the 
negative sign problem. 
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FIG. 5: (Color online) LTRG and TMRG results of the quantum XY 
spin chain, (a) Relative error of the free energy per site Sf. (b) The 
energy per site e. The inset shows the variation of (e - e )/e with 
inverse temperature ft for various D c . 
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FIG. 6: (Color online) (a) Specific heat as a function of temperature 
(T = 1 IP) of the quantum XY spin chain. The inset shows the low 
temperature results for D c = 100 and 200, along with the TMRG 
data (M = 200) for a comparison, (b) Energy per site of the 2D spin- 
1/2 Heisenberg antiferromagnet on a honeycomb lattice for different 
staggered magnetic fields. The QMC results are obtained by using 
the ALPS library (Ut 

ture (T/J =* 0.008). The TMRG results with states M = 200 
are also included, showing that both numerical methods have 
the comparable accuracy. 

To examine the scalability of the LTRG algorithm, we also 
calculate the energy per site of a 2D spin- 1/2 Heisenberg an- 
tiferromagnetic model on a honeycomb lattice, whose Hamil- 
tonian is H = JZ<ij>^i ■ + K^-lfS*, where (-1) 1 '' 1 
denotes the parity of the lattice and h s is a staggered mag- 
netic field, as shown in Fig. |6](b). A pronounced agreement 
between LTRG and quantum Monte Carlo (QMC) results is 
clearly seen. 

In summary, we have proposed a linearized TRG algorithm 
to calculate the thermodynamic properties of low dimensional 
quantum lattice models, and obtained very accurate results. 
The LTRG algorithm can be readily generalized to fermion 
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